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Weakly focusing transparent media provide possible means for guided 
transmission of coherent light beams with relatively small loss. The scalar 
wave equation for the eigenmodes of propagation in such a medium is 
formally identical with Schroedinger' s wave equation. Hence, the methods 
used in the solution of quantum-mechanical problems, such as the Wentzel- 
Kramers-Brillouin (WKB) approximation, are immediately applicable to 
this problem. Solutions for the eigenmodes and eigenvalues in the case of 
focusing in one dimension are given, and the Poschl-Teller medium, whose 
index varies as 

n = n [l - (a/2) tan 2 v x] 

is discussed in some detail. In addition, the relationship between the wave 
solutions and geometrical (ray) optics is examined. 

I. INTRODUCTION 

Weakly focusing transparent media, exemplified by the gas lens, 1 pro- 
vide a possible means for guided transmission of coherent light beams 
with relatively small loss. The optics of a medium whose refractive index 
decreases quadratically away from some spatial axis (the z axis, say) has 
been the subject of much discussion in the literature. 2 In this paper, we 
would like to examine the properties of more general guiding media. Some 
work on this problem has been carried out recently by S. E. Miller. 3 

Consider a light beam traveling paraxially in the z direction in space, 
guided in a weakly focusing transparent medium. For simplicity, we 
shall consider focusing only in one of the two transverse dimensions. 
Thus, the refractive index of the medium and the electromagnetic field 
will be assumed to be dependent only on the x dimension, and to be 
independent of y. 

Our concern is with transparent media whose index of refraction has 
the form 

n = ?i„(l - %f(x)). 
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We make the following assumptions about /(.r) 

f(x) « 1 in the range of interest (1) 

U a ° (2) 

dx 2 

/(0)=0; f(x)=f(-x). (3) 

We will investigate both geometrical and physical optics of such a 
medium in the realm of validity of the paraxial ray equation. Of the 
assumptions, (1) implies only gradual changes of index, (2) insures 
focusing properties, and is used in the approximation procedure of 
Section Vll, while (3) is made only for mathematical convenience and 
can easily be relaxed. 

II. GEOMETRICAL (RAY) OPTICS 

The well-known paraxial ray equation 

(tx = Idn = _ldf 
dz 2 n dx 2 dx 



(4) 



has the following general solution. 
If we let 



dx 

v= Tz> 



(5) 



i.e., p is the slope of the ray path; then (4) becomes 

dp = _ldf 
P dx 2dx' 

Hence, 

p 2 _|_ f = cons t = £. (6) 

Inserting (6) in (5), and solving, we get 

2 = / /* ti \ + const - ( ? ) 

Note from (6), (5) and (3) that ±£* is the slope of the ray path as it 
crosses the axis x = 0. 

III. PHYSICAL OPTICS 

In the paraxial ray approximation, the electromagnetic fields are 
always very nearly perpendicular to the "direction" of propagation, 
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which we take to be the z axis. Hence, we can use the scalar wave equa- 
tion. For harmonic [exp (—itri)] fields that are independent of y, we 
arrive at the wave equation 

^ + P 2 + fcV(l - /(*))* = [* - «/c] (8) 

dx 2 dz 2 

where we have used (1) to expand n 2 . We look for solutions to (8) with a 
propagation constant /3 in the z direction, i.e., 

* a '""■ % - -** (9) 



yielding, for (8) 



p 2 + [/3„ 2 - /3 2 - tff(x))+ = (10) 

dx 2 



where we have substituted n 2 k 2 = O 2 = (2?r/X) 2 . With the further 
substitution 



i - 2 /Po 2 = a (id 

(10) becomes 

^ + ft 2 (£ -/(*))* = 0. (12) 

Equation (12) is Schroedinger's wave equation for a particle in a one 
dimensional potential well. It has a sequence of eigensolutions yp m and 
corresponding eigenvales £ m . The eigenfunctions \j/ m here represent the 
transverse distributions of the propagating field modes, while the 
eigenvalues £ m give the propagation constants. In accord with (1), 
the eigenvalues will be much smaller than unity, and so we can find 
the propagation constants /3 m from (11) as 



(i-tK 



(13) 



If we use the notation of the classical ray path equation (G), and set 

V = P(x) = VZ - f(x) (14) 

then (12) becomes 

g + A.W-0. (15) 

We note from (2) and (3) that f(x) increases monotonically with 
I x I. Consider the points x = ±A, where 
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/Oc)-/(A)-«. (16) 

For | re | < A, p 2 is positive, and \p has oscillatory behavior. For \x\ > A, 
p 2 is negative, and yp has decreasing exponential behavior. The points 
| x | = ±A are inflection points of \p. In the ray solution, the points 
j re | = A, where p = 0, are the turning points, where the ray has zero 
slope and, correspondingly, maximum excursion from the axis. 

IV. THE WKB APPROXIMATION 

A well-known solution to (15), valid in the range | x | < A if p can 
be considered approximately constant over a few cycles of the oscilla- 
tory behavior of \f/, is 



\f/ cc p * cos <ft, / pdx + a>. (17) 

With the symmetry of / assumed in (3), one can show directly from (12) 
that \f/ must be either symmetrical or antisymmetrical in x. Thus, a 
must be a multiple of tt/2. The eigenvalues of £ are then determined by 
the condition that \p must be matched through the turning point x = A 
(where another approximate solution to (15) is necessary) into the 
decreasing exponential solution for x > A. Discussions of this problem 
may be found in most quantum mechanics texts. 5 If f(x) can be assumed 
linear in x over a suitable region near A , then the asymptotic formulae 
for connection through x = A are 



\j/ cc p ' cos ( 6 — - j 


x < A 


(18a) 


-*(|*)VtW +/-!(•)] 


x < A 


(18b) 


-*(!^y [-/»(») +/-iw] 


x > A 


(18c) 



where 



ip -} exp(-0) x>A (18d) 



9 = e(x) = I ft, f \p(x)\dx'\. (18e) 

J A 



The phase of ^ as it approaches the turning point must be as in (18a). 
Otherwise the connection through the turning point would give rise to 
the increasing exponential in the region x > A, and this is unallowable. 
Expressions (18a) and (17) for ^ can be equated only when 
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ft, f Pdx= (TO+|)g (19) 

Jo Z 

where m is an integer. Even and odd symmetry solutions correspond 
respectively to even or odd values of m. Equation (19) is Bohr's quanti- 
zation rule, and its solutions give the eigenvalues £ m . The eigenfunctions 
are then given by (18). 

V. RELATION BETWEEN WKB AND GEOMETRICAL RAY PATHS 

Referring to (19), let 

Kk) =— f Vdx. (20) 

7T •'O 

For successive eigensolutions, I changes by unity. Now, 

S-fiff* (21) 

where we have used (14); i.e., p = V£ — f(x). For a change of / by 
unity, the propagation constant (13) changes approximately by 

8 *~ d i dt/jr; r^/p (22) 

J o 

Now from the ray solution (7 ) , the z distance a ray with the same value 
of £ takes to go through one complete cycle of its transverse oscillatory 
motion is 



z_ = 4 / dx/p (23) 

Ja 



•A 
'0 

hence, putting this in (22), we find 

(80)z„ = 2tt. (24) 

Hence, adjacent modes undergo a relative phase shift of approximately 
2x in the same distance the corresponding geometrical ray takes to 
complete one transverse cycle. 

We can make the same argument more general by examining the be- 
havior of solutions of (15) which are sums of component solutions (17) 
with nearby values of f. We look for trajectories x(z) such that the 
phase differences between components remain constant. Then if we make 
up a wave packet, the packet will follow one of these trajectories. 

Putting the z dependence back into (17), and keeping, for example, 
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only the positive imaginary part of the cosine term, we have 

if, cc p _i jexp i I 02 + ft, j pdx + a\\. (25) 

The phase we have to maintain constant over a range of £ is thus 

0z+ P„ [ pdx + a. (26) 

Jo 

Differentiating with respect to £, setting the result equal to zero, and 
using (13) and (14), we have 

*=/*- + «-£• (27) 

Jo p ft. d£ 

To complete the picture, we note that in order to make a wave packet 
initially concentrated at some point (x„ , z ) should take all the phases 
initially equal at that point, so that 



a = -$z - ft, / p 
Jo 



2 da f x ° dx 



dx + const 



Po d£ Jo p 

Hence, the wave packet trajectories of stationary phase are given by 

z - Zo = r*t. (28) 

J x P 

This is identical with the ray path (7). In this last discussion we have 
neglected the discrete nature of the eigenvalues. Such an approximation 
should be valid for reasonably large mode numbers m. 

VI. THE POSCHL-TELLER POTENTIAL 

A number of functions f(x) yield analytically integrable equations. 
The square law medium f(x) = ix/bf and the square well medium 
f{x) — for | x | < A, f(x) — > oo for x > A, suitably joined, are media 
with well-known solutions both for the geometrical and physical optics 
equations. An interesting function, which in a slightly more general 
form goes by the name of the Poschl-Teller potential 6 in quantum 
mechanics, is the function 

f(x) = a tan 2 (nx); -t/2jj < X < r/2 V . (29) 

Near x = 
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f(x) = arfx 2 (30) 

and there are impenetrable boundaries at x — ±7t/2tj. For this function 
the ray equation (7), the wave equation (12), and the WKB approxi- 
mation (17) can be directly integrated and compared. As the parameters 
a and t\ are changed, the function varies smoothly between a square 
well and a square law type. 
Consider first the ray equation, 

/dx 
i - a tan' „]• + C0DSt (31) 

This integrates to [Burrington % 258] 

z =—4—77 sin" 1 ([1 + «/£]* sin t\x) + const. (32) 

If we set the constant equal to zero, then 2 = when x = 0. This expres- 
sion generates a ray path, with x taking on values between the turning 
points 

x= ±A = irf ' sin -1 [1 + a/*] - * (33) 

while the period of oscillation is 

z~ - 2*/nte + a]K (34) 

Note that at the turning points, f(x) = £. We see that for a » £, i)x 
is always small, and (32) reduces to 

x = ($/tf 2 tt) J sin (ijaa) 

which is appropriate to a square law medium, while for a « £, the ray 
travels back and forth with constant slope £* between reflecting walls 
separated by ir/r]. 

Poschl and Teller have found exact solutions for the wave equation 
(12) with this potential. They have shown that the eigenvalues follow 
the simple law 

£ m = (v 2 /0o) (™ + 2ma + a) (35a) 

where a is the positive root of the equation 

a(a- 1) = atf/rj 2 . (35b) 

The corresponding eigenfunctions are 

even 
m odd 

\p m = cos" (rji-) ^ c k sin* (tjz) (36) 

A=0,even 
Ik— l.odd 
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where 

c/t+2 _ (ra — k)(m + k + 2a) 
c k (k + l)(k + 2) 

From (35) and (13), we find that the propagation constants of the 
modes follow the law 



ft» = & - (va/2Po) U + 2m + ^V 



(37) 



To examine the square law limit, we set a » £. Then [see (35)], a « 
a%/ V » 1, 

and, from (37) 

j8««0 o - r,aHm + h) (38) 

while from (36) 

* m PX exp (-±r,<x%T 2 )H m ( V a%x 2 ). (39) 

To examine the square well limit, we set a « £. Then a tt 1, and 

/3 m -+ft-M'?7/3 )(m+l) 2 (40) 

^ m — > cos [(ra + l)i?x] for ra even 
!/-„, — » sin [(ra + 1)ij.t] for m odd. 



(41) 



Returning now to the general solution for the propagation constants, 
note from (37 ) that the average of the propagation constant differences 
between the rath mode and its two neighbors is given by 

5/3 = \ 03 m _ a - m+1 ) = (v/Po)(m + a). (42) 

If we look at the ray equation, and insert the values of a and £ m (35) 
into (34), we find 

z„ = 2ir/8j3. (43) 

This is a more precise version of (24). 

Using the Poschl-Teller potential, the WKB result (19) for the 
eigenvalues £ ro can also be integrated. Referring to (21), (23), and (34), 
we have 

~ = (A/t) T dx/ V = /W4tt 
a\ Jo 
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Hence, 

/ = &(£ + «)*A - ft«V* (44) 

Here the constant of integration is determined by the condition that 
I = at £ = 0. Thus, (19) [see also (20)] is 

fttt- + «)Vl = ^Vu + (m + |) (45) 

or 

6-- Kf/fcX* +*) +«*f-a 

(46) 

= (v/Po) 2 [(m + £) 2 + (2m + 1) Va(o - 1)] 

where we have used (35b). We see that the WKB result (46) agrees 
with the exact result (35a) in the limit of the square law medium, 
a » 1. In the square well limit, a — » 1, the WKB method gives eigen- 
values proportional to (m + \) 2 rather than (ra + 1) as does the exact 
result (35). This last discrepancy can be traced to incorrect matching 
of the boundary conditions by the WKB method, since here the WKB 
wave function (17) is the correct one. The factor \ in ra + \ arises from 
matching the boundary condition, and in fact, we can get the WKB 
answer to equal the exact one for all a and ra here by the artifice of 
replacing m + £ in (46) by 

m + a - \/a(a - 1) = m + a - otpo/v- (47) 

The number added to m varies between the value \ for large a to one 
as a approaches unity. 

VII. FURTHER APPLICATION OF THE WKB METHOD 

The results of Section VI give us reasonable faith in the WKB method 
of obtaining eigenvalues even at small values of ra. The number we add 
to ra may however, take on values between \ and 1 depending in general 
on m as well as on the form of the potential. 

Let us consider the WKB result further. Equation (19), which we 
repeat here for convenience, may be written 

(Jo f m Vtm - Six) dx = (m + i) ^ (48) 

where 

f(A m ) = £„, . 
We can see two simple bounds. For functions that satisfy (2) and (3), 
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we have that in the range ^ x ^ A m , 

U^Zm- /(*) ^ 6. (l - jj (49) 

hence, performing the integral for the two bounds, (48) gives 

(3 VbnA m ^ (to + h) | ^ |/3. VMm- 

We can write this as 

VT m A m = r(m + |)ir/2A, (50) 

where 

1 ^ r ^ 1.5. 
Note that 

r~ l = ^ jT V* - /(*) ifa. (51 ) 

If we further remember that we should replace the % in to + £ by some 
number between £ and 1 which we can denote by s, we have 

VU™ - Km + sV2& (52) 

with 

lgr^ 1.5 and 0.5 ^ s ^ 1. 

Let us examine the implications of this formula in a definite example. 
We take the potential to be a pure power law 



fix) - ( I ) . (53) 



We need the value of A m ; 



f(A m ) = £m therefore, A m = bf« ! *. 
Then (52) is easily solved, yielding 

In this case, examination of (51 ) reveals that r is given by 

r _1 = f Vl - y** dy (55) 

Ja 
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so that r is independent of m, and varies between 

4 

for n = 1 and 1 for n — * °o . 

7T 

Further, we know that s = £ for n = 1 (square law), and s = 1 for 
n — » oo (square well). As a guess, we would be tempted to try 

n 

s = 



n + 1 
as a suitable interpolation. Then we have as our final result 

hence for the phase constant [repeating (13)] 

fi m = fc(l - f m /2). (57) 
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